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ABSTRACT 

We discuss the systematic effects arising from tlie cosmological redshift- space (geo- 
metric) distortion on the statistical analysis of isodensity contour using high-redshift 
catalogs. Especially, we present a simple theoretical model for isodensity statistics in 
cosmological redshift-space, as a generalization of the work by Matsubara (1996). The 
statistical quantities considered here are the two- and three-dimensional genus of iso- 
density contour, the area of surface, the length of contour intersecting with a plane and 
the number of the crossing points of isodensity contour on a line. We give useful analytic 
formulae for the isodensity statistics, which take into account the corrections from the 
geometric distortion, the nonlinear clustering and the nonlinear velocity distortion phe- 
nomenologically. We then demonstrate how the geometric distortion and the nonlinear 
corrections alter shapes of the statistical quantities on the basis of plausible cosmological 
models. Our results show that the nonlinear correction can be sensitive to a choice of 
the redshift-space coordinate as increasing the redshift. The low-dimensional quantities 
such as two-dimensional genus systematically yield anisotropy due to the geometric and 
velocity distortions and their angle-dependence shows the 10 ~ 20% difference of ampli- 
tude. Sensitivity for typical high-redshift samples are estimated in an analytic manner, 
and the influence of the light-cone effect for the isodensity statistics is also discussed. A 
simple estimation suggests that the systematic effects of geometric and redshift-space 
distortions can become comparable or could be dominant to the statistical error of deep 
cluster samples and future high-redshift galaxy surveys. These systematic effects might 
be a useful tool in probing the cosmological model of our universe. 

Subject headings: cosmology: theory - dark matter - large-scale structure of universe - 
galaxies : distances and redshifts - methods: statistical 
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1. Introduction 

Recent observations of high-redshift objects have revealed dynamical aspects of the evolving 
universe. Measurements of Lyman-break galaxies provide evidences for the early galaxy assembly 
at the epoch z ~ 3, which put a serious constraint on the structure formation scenario (Steidal et 
al. 1998; Giavalisco et al. 1998). Clusters up to z ~ 1 have been observed, and the abundance 
suggests the universe with a low density parameter 0, ~ 0.3 (c.g, Carlbcrg et al. 1997). Moreover, 
the upcoming wide-field surveys such as the Sloan Digital Sky Survey (SDSS) and the Two-Degree 
Field (2dF) survey promise to extend the observable scale of the universe and these enormous 
data will precisely reveal the clustering feature of the high-redshift quasars and galaxies. Since 
the analysis combining different redshift data sets has a potential to break the degeneracy of the 
cosmological parameters, as well as to provide the strong constraints on the structure formation 
scenario, statistical studies using the high-redshift objects arc of great importance as a complemen- 
tary to the measurements of cosmic microwave background anisotropy and the detection of cosmic 
shear. 

Recently several authors have discussed various cosmological effects arising from proper obser- 
vations of the high-redshift objects (e.g, Suto et al. 1999 and references therein). Apart from the 
biasing of the luminous objects, the main effects are summarized as follows: 

1. geometric distortion : no one can obtain the correct three-dimensional map of the high-redshift 
objects without knowing the correct cosmological model of the universe. The uncertainty of 
the cosmological model may affect proper evaluation of the clustering pattern in the cosmo- 
logical redshift space (Alcock Sz Paczyhski 1979). 

2. redshift-space distortion : peculiar motion of a cosmological object affects the estimation 

of the distance in a redshift survey, which causes the linear redshift-space distortion due to 
the bulk motion and the nonlinear velocity distortion due to the virializcd random motion 
(finger-of-God effect) (Davis & Peebles 1983; Kaiser 1987; Suto &: Suginohara 1991; Hamilton 
1998 for a review). 

3. light-cone effect : because an entire sample does not consist of objects on a constant-time 
hypersurface but on a light cone, the redshift evolution of the objects inevitably contaminates 
the data. Therefore, in a wide and deep survey, it becomes important to incorporate the 
light-cone effect into statistical quantities, for example, two-point correlation function and 
the power spectrum (Matarrese et al. 1997; Yamamoto & Suto 1999). 

For the popular two-point statistics, the significance of these effects have been studied in 
detail and the robustness of the theoretical prediction is checked by the numerical simulation 
(Ballinger et al. 1996; Matsubara & Suto 1996; Magira et al. 2000; Hamana et al. 2000). Using 
the theoretical formulae, the feasibility of cosmological test with geometric distortion has been 
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investigated (Yamamoto, Nishioka & Taruya 2000). Further, the analysis is developed to the 
higher-order moments on a light-cone (Matsubara et al. 1997). 

Here we extend these analyses to the isodensity statistics. That is, we investigate the influences 
of the above cosmological effects on the isodensity statistics, which have not been ever clarified. The 
isodensity statistics Uke the genus G3 and area statistic N-^ have progressively become important 
as a quantitative measure of the topology of the large-scale structure. Analytical expressions of the 
isodensity statistics are derived for the Gaussian random field and some cases of non-Gaussian field 
in real space (e.g, Bardeen et al. 1986; Hamilton et al. 1986; Matsubara 1994, 2000; Matsubara & 
Yokoyama 1996). With these analytic formulae, the isodensity statistics computed from the redshift 
surveys of galaxies and/or the simulated dark matter catalogues are comprehensively understood 
(Colley et al. 2000 and references therein). 

While most authors of these works have focused on the low-redshift objects except for the 
cosmic microwave background anisotropy, the analysis using the high-rcdshift samples is expected 
to provide an important clue on the origin and evolution of proto-galaxies and proto-clusters. In 
fact. Park et al. (2000) have recently analyzed the two-dimensional genus using the galaxies taken 
from the Hubble Deep Field and quantified the feature of early galaxy assembly. Remarkably, their 
samples are widely distributed over the redshift range ^ 2; ^ 3, which revealed an unprecedentedly 
deep view of galaxy clustering. Although the analysis is still limited to the projected density field 
because of the large error of photometric redshift, precise measurement of redshift will enable us 
to quantify the various isodensity statistics defined in the three-dimensional space, as well as the 
quantities in the one- and two-dimensional space in near future. 

Practically, the observed distribution of high-redshift objects suffers from the cosmological 
effects, which should be correctly taken into account when comparing such data with theoretical 
prediction. Then naive but important questions are as follows: how the cosmological effects affect 
the statistical study of isodensity contour; whether these effects are observed as signal or noise. 
Purpose of this paper is to address these issues by using a simple theoretical formulae for the 
isodensity statistics, in which the cosmological effects are faithfully incorporated. Partially, the 
influence of the linear (redshift-space) distortion is investigated by Matsubara (1996), though the 
effect turns out to be small. In his paper the corrections from the nonlinear clustering and the 
nonlinear velocity distortion, which become influential on small length scales, are not investigated, 
and neither is the geometric distortion. 

In this paper, we generalize the Matsubara's results to those incorporating the geometric 
distortion effect and the nonlinear correction taking account of the smoothing effect. With the ex- 
tended formulae, we examine how the geometric and redshift-space distortions affect the isodensity 
statistics qualitatively and quantitatively. Wc will show that the nonlinear effect can be important 
even for a case with large smoothing scale because the scaling effect of the geometric distortion 
causes apparent enhancement of the nonlinear effects, which depends on choice of the cosmological 
redshift space. In addition we also discuss the sampling noise and the effect of finite redshift depth 
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for typical high-redshift samples by using our extended formulae. 

We organized the paper as follows. In section 2, the cosmological redshift-space coordinate 
is introduced, and the statistical quantities incorporating the nonlinear corrections are described. 
Then we consider the isodcnsity statistics in cosmological redshift space in section 3. Based on the 
analytic formulae for the Gaussian field in appendix A, we present a nonlinear model of isodcnsity 
statistics. After demonstrating predictions of the model, sensitivity of the isodcnsity statistics and 
the effect of the light-cone are discussed in section 4. The final section 5 is devoted to summary 
and conclusion. 



2. Density field in cosmological redshift space 
2.1. Redshift-space coordinate 

We start with the standard form of Friedman-Robertson- Walker metric 

ds'^ = -c^dt^ + a'^it) [dr^ + fK{rf {dO^ + sin^ edcf)'^)] , (1) 
where Jk{t) is the comoving angular diameter distance given by 



fK{r) = < 
with the spatial curvature K 



sin {^/Kr)/^^K (K > 0) 

r (i^ = 0) , (2) 

sinh(x/^r)/x/^ (iC < 0) 



K= (^Vcj^o + Ao-l), (3) 



{z) = ^ ; H{z) = iJo V^^o(l + zr + (1 - - Ao)(l + z)^ + Aq. (4) 

Jo J^\z ) 



and the comoving radial distance 

cdz' 

'0 



The quantities, VLq^ Aq, and Hq^ respectively, denote the density parameter, cosmological constant, 
and the Hubble parameter, at the present epoch z = 0. 

Because we are interested in a redshift survey, suppose that we could only know radial distances 
of objects in terms of the redshift z and separations between pairs of objects in terms of the 
angular diameter distances. This means that the relation between the comoving real space and the 
observable redshift space inherently depends on an assumption of the cosmological model of the 
universe, and that the radial component of peculiar velocities of the objects apparently alters the 
distribution in the observed redshift space. 



Hereafter, unless otherwise stated, we will assume that the sample of objects located in the 
range z — dz and z + dz with dz <^ z whose angular separations are of the order of dz. Then, the 
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distant-observer approximation is entirely valid. Within this approximation, we consider the object 

located at r = (ri,r2,r3) in the comoving real space. For dcfinitcness, we choose the third-axis 
as the hnc-of-sight direction. Denoting the corresponding coordinate in the observed cosmological 
redshift space by s = (si, S2, S3), the relation between s and r is given by 



Si = 



<7i 


_ ri 




c±iz) 


Q2 


_ r-a 




c±{z) 




1 




c\\iz) 



c±{z) Cj_[z) 

where is the linc-of-sight component of the peculiar velocity field. The observed redshift space 
s is obtained from the new coordinate system q = {qi,q2,Q3) divided by the dimcnsionlcss factors 
Cx and C||. Since we must assume a distance-redshift relation s = s{z) to plot a map, the factors 
Cx and C|| depend on the assumption of the cosmological redshift space. Frequently, one adopts 
s = ■SHb('Z)) the radial-distance relation extrapolating the Hubble law, and/or s = s-E^iz), the 
relation in Einstein-de Sitter universe: 

C Z 

SHbiz) = — , (6) 



^0 

c 



s..iz) = l^fl-—^). (7) 



Then, the geometric factors become 



c 



for the case s = SnbC^) (Matsubara &; Suto 1996), and 
for the case s = SEdiz) (Balhnger et al. 1996), respectively. 



2.2. Field correlations 

Recall that the number density fluctuation of luminous objects is observed in the cosmological 
redshift space s. We can evaluate the fluctuation (5(*)(s) smoothed over the spherical radius Rs : 

5^'^ (s; Rs) = J (fs' W{s - s'; Rs)5^'^ (s') 

= / ^W{ksRs)5^'HK)e-''^'-^ (10) 
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where the quantities W and d^^") are the Fourier transform of the window function W and the 
density field 5^^^ , respectively, is the Fourier vector corresponding to the s-space. In this paper, 
we adopt the Gaussian window function, i.e., W{y) = exp(— The density fluctuations defined 
in q- and r-space are merely related to 5^^') through the transformation (5). 

For our purpose, it is convenient to introduce the new field variables: 

n, = di6^'\ Cij = dAS^'^ ; = 1,2, 3), (11) 

where di is the derivative with respect to the cosmological redshift-space coordinate Sj. Since 
correlation functions in the redshift space are inherently anisotropic for the special direction of 
the line-of-sight, the variance and the cross- moments of the variables (11) have the following non- 
vanishing quantities: 

= al (12) 

2 

iiVs?) = al, , in.nj) = ^d,j = - {Oj6^^^), (13) 

2 

), (14) 

where (• • •) denotes ensemble average over the smoothed density fields. Subscripts / and j run from 
1 to 2. Tlic quantities ctq, cJi h, o"i_j^ and o"2,x are related to the power spectrum p'-"'^')^ defined in 
the cosmological redshift space s. Substitution of (10) into (12), (13) and (14) yields 

W'ikM p(-'^>(fe„|i,fc„x;^), (15) 
<^l,{Rs) = l^w'{ksRs)P^-''\K^uks,±;z)kl^,, (16) 

(27 

with 



<ylARs) = J ^w'iksRs) P'-"''\K,n,hx,z) (18) 



ks — \/ k'^u + k^ 



where ks^^^{ks^±) denotes the wave number parallel(perpendicular) to the line-of-sight direction. 

We notice that the power spectrum P^"'^^(A;3 n, A;j,_x; -2) is written by the power spectrum defined 
in the q-space, P^^\k\\, kj_;z), where, and k^ arc respectively the line-of-sight and its orthogonal 
component of the Fourier vector corresponding to the q-space. Using the facts that fc^ji = C||fc|| and 
^s,x = Cj_kj_, we obtain 

p<-)(fc.,,„V) = -V^^^^(7^>^)- (19) 

Cm Cx V C|| Cj_ J 
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Thus, all of the field correlations can be expressed in terms of the quantities defined in the q-space: 

W\kR,s) P^'\h,k^;z), (20) 
W\kR,s)P^'\k„k^;z)kl (21) 







J (27r)3 


(^liiiRs) 


= J 


<^lARs) 


= cliz) j 


<^lARs) 


= ci{z)j 



(27r)3 

^ W\kR,^) P^'\k,M\z) kl, (22) 
W\kR,s) P^'\h,k^;z) ki, (23) 



^ (27r)3 

where k = ■^kf^ + k'^. The quantity R^s is related to the smoothing radius R. 



ReS = c±Rs 



fell \ L f 



(24) 



Note that the spherically symmetric smoothing in the s-space corresponds to the anisotropic filtering 
in the q-space in general cases with ^ c±. This can be regarded as an effect of the geometric 
distortion on the one-point statistics, as well as the isodensity statistics. 

Now we need a model for the power spectrum, P(*)(/C||, kj^; z) (or p(<=''d))^ including the rcdshift- 
space distortion and nonlinear clustering. Phenomenologically, but fairly accurate fitting form of 
the power spectrum has been recently proposed by several authors (Cole et al. 1995; Ballinger 
et al. 1996; Suto et al. 1999), and its validity is successfully confirmed by N-body simulation 
(Magira et al. 2000). Hereafter, we assume that the observed number density fluctuation is simply 
proportional to the mass density contrast by introducing the biasing factor 6 as the function of 
redshift 

d^^Hr,z) = b{z)5(^Ur,z), (25) 
in real space. Then the power spectrum of the observed density fluctuations is written as 

Pif{k„k^;z) = b\z)Pj^Uk) [l + m^^'Y D[kl^^Az)] , (26) 

where Pnilssik) is the power spectrum of mass fluctuation in real space, cjp denotes the peculiar 
velocity dispersion, and we defined fi = k^Jk. The factor [1 + describes the linear theory of 

redshift-space distortion (Kaiser 1987) and (3 is obtained through the differentiation of the linear 
growth factor Dj^{z): 



b{z) din a b{z) [ ^ ' 70 
with 

-"0 Wi , n ■ — I tin 



(27) 



^^'^=^H{z)) (1 + ^)^°' ^(^)=(,^) (28) 
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The last term in (26) is the damping factor due to the finger-of-God effect: 

which corresponds to the assumption of the exponential form for the distribution function of the 
pair-wise peculiar velocity in real space. To be specific, in this paper, we adopt the fitting formulae 
of Peacock & Dodds (1996) and Mo et al. (1997) for the nonlinear mass power spectrum Pmass{k) 
and the peculiar velocity dispersion Cp, respectively. Then the field correlations, a^, (Ti_||, cti^x and 
(72,_L can be computed through the definitions (20-23). 



3. Isodensity statistics in cosmological redshift space 

In this section, we consider simple theoretical models for isodensity statistics in cosmological 
redshift space. After describing definitions of the isodensity statistics, plausible nonlinear formulae 
are presented, based on the Gaussian results in appendix A. Influences of the nonlinear corrections 
and the geometric distortion are quantitatively examined in familiar cold dark matter cosmological 
models. 



3.1. Definitions 

First consider the three-dimensional genus curve. The genus G3 is defined by —1/2 times the 
Euler characteristics of two-dimensional surfaces. The Euler number density is given by (number of 
holes) — (number of isolated regions). Equivalently, this can be expressed by (number of maxima) 
-I- (number of minima) — (number of saddle points) of the isodensity contour with respect to 
some fixed direction irrespective of the isotropy of the random field. At a given density threshold 

(5(^)/(7o, the genus G3 is defined by (Bardeen et al. 1986) 

G^i\^) = -\ {Sd{S^''^ - i^cj,)5D{m)5D{ri.) 1%! (CnG^ - C12)) , (30) 

where 5d{-) denotes the Dirac's delta-function. 

Genus is also defined for a two-dimensional random field. Consider the field in a two-dimensional 
fiat plane S and its isodensity contours. The two-dimensional genus G2 can be computed by count- 
ing the number of the isolated high(low)-dcnsity regions surrounded by the isodensity contours 
(Melott et al. 1989). Notice that the quantity G2 depends on a configuration of the plane S be- 
cause the distortion induced by the geometric and peculiar velocity effects makes the random field 
anisotropy (Matsubara 1996). Denoting the angle between the line-of-sight direction and the plane 
S hy 6 (see Fig.l), the two-dimensional genus per unit area, G2, is 

g'^\u,9) = ~ (<5d(<5(^) - uao)SD{Vi) |r/2 sin^ + r?3 cos0| Cn) • (31) 
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Next consider the area statistic N^. This was proposed by Ryden (1988) and Ryden et al. 

(1989) as a complementary to the genus statistics. Given a density height u, the quantity A^3 
quantifies the mean surface area of the isodensity contours per unit volume : 

= {Sd{S^'^ - i^cTo)VW+W+vl ) ■ (32) 

For a homogeneous and isotropic random field, A'^^*''(i^) equals A/tt times the mean length of the 
isodensity contour with same height u per unit area on a flat plane and twice the mean number of 
the crossing points per unit length on a straight line, which we denote by N2 and A^i, respectively 
(Ryden 1988). In cosmological redshift space, these statistics obviously depend on the angle 9 
between the line-of-sight direction and the plane S or the line L (Fig.l). Therefore we need treat 
N2 and A^i, separately, which can be defined as 

N^'\iy,9) = (sd (S^'^ - i^To) \lri + iv2 sin ^ + rys cos ^)2 ^ , (33) 
Ni'\iy,0) = (^SoiS^'^ -iy<7o) \V2 sine + ris cos 9\y (34) 



3.2. Modeling the nonlinear density field and isodensity statistics 

Provided the distribution function for the field variables S^^\ rii and Qj, the expressions for the 
isodensity statistics defined in previous subsection can be derived analytically. For the homogeneous 
but anisotropic Gaussian random fields, we easily obtain the simple analytic formulae (appendix 
A). The influence of geometric distortion is naturally incorporated into them. 

The Gaussian results, however, are rather idealistic and cannot be used beyond the linear 
regime. Properly, we must take into account non-Gaussianity of the density distribution function 
developed by the gravitational clustering in the cosmological redshift space. Since no rigorous treat- 
ment of nonlinear evolution has been found, we need first adopt a reliable model of the distribution 
function V^S^^^ , rji, Qj). A class of models for the nonlinear density field have been discussed, whose 
statistical property is characterized by a monotonic function F (e.g. Coles & Jones 1991): 

<5W(s)=F[</.(s)], (35) 

where the scalar field (f){s) is a Gaussian random variable with zero mean and the unit variance. In 
this case, it is easy to derive analytic expressions of isodensity statistics from the Gaussian results 
in appendix A. This is partially done by Matsubara &; Yokoyama (1996) for the genus curve G3. 
We generalize their result to the other statistics G2, N^, N2 and Ni. In appendix C, assuming the 
functional form (35), the analytic formulae for isodensity statistics in the non-Gaussian case are 
presented, which are simply obtained by transformation of the Gaussian results. 

The plausible functional form (35) for the nonlinear density distribution has been investigated 
by several numerical works. They showed that the probability density distribution of one-point 
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statistics can be empirically fitted by the log-normal distribution function to a good accuracy 
(Kofman et al. 1994; Bernardeau Sz Kofman 1995; Taylor & Watts 2000). This is also the case 
in the redshift space (Hui et al. 2000). Observationally, the log-normal distribution function has 
been used as a model of the galaxy distribution ever since Hubble (1934) (e.g, Hamilton 1985; 
Bouchet et al. 1993; Kofman et al. 1994). Hence, we might use the log-normal function as a good 
approximation for the actual observed distribution. 



The log-normal distribution is characterized by 

1 



m = 



exp [ai,ti(f>\ - 1, 



where 



Then we obtain the one-point distribution function 



exp 



{ln(l + <5W)+c7L/2} 



21 



_d6^ 
1+^' 



and 



ln(l + uao) + 



(36) 

(37) 

(38) 
(39) 



Using (39) and (C2)-(C6), the analytic expressions for log-normal distribution are obtained as 
follows : 



47r2(l + (72)3/2 a 



3 

LN 



1 



X exp 



ln(l + t/ao) + ag^/2 |'' 
[ln(l + i^a,) + 



(40) 



for the three-dimensional genus, and 



ln(l + i/cTo) -f- 



X exp < — 



[ln(l + z/ao) + ag,/2]2 

2(72 



for the two-dimensional genus. The quantities, N3, N2, and A'"i reduce to 



1 



0-1 



2v^7r(l + C72) (Tln 

1 cri,± 
2v^7r(l + (72) (Tln 

2y7r(l + (72) (Tln 



exp < - 



exp < - 



exp 



[ln( l + z/(7o) + al^/2f 
2fT2^ 

[ln( l + lyg,) + c7g^/2]^ 
2ct2^ 

[ln(l + lygp) + (7^/2]- 
2ct2, 



i^(7,^), 
i^(7,^), 



(41) 

(42) 
(43) 
(44) 



where the functions J, K and L are respectively given by (Bl), (B2) and (B4) in Appendix B, and 
we introduced the distortion parameter 7 given by 

^ = ^2^. (45) 

Notice that the parameter 7 equals unity when C|[ = = 1 and the power spectrum is isotropic 
P^^\k\\,kj_) = P^^\k) from equations (21) and (22). When the geometric distortion is effective, the 
factors C|| and are not equal to unity. Further, when the linear distortion or the finger-of-God 
effect is influential, the power spectrum is no longer isotropic. Thus the deviation of 7 from unity 
quantifies the collected distortion effects, and represents anisotropics of the isodensity contours. 
As will be described in section 3.3, 7 is affected by the three distortion effects: the geometric 
distortion, the linear distortion, and the finger-of-God effect. The linear distortion increases 7. On 
the contrary the finger-of-God effect decreases 7. As for the geometric distortion, this depends on 
the choice of cosmological redshift space s{z) (see cqs.[21] and [22]). When we adopt s = SHb(^)) 
7 is decreased due to the geometric distortion. Conversely when adopting s = s^^, 7 is increased. 
These features of the geometric distortion arc based on the assumption that the real universe is 
the Lambda CDM model (see below). The strength of each distortion effect depends on smoothing 
scale and redshift, which yields the scale-/time-dependence of the value 7. 

The above formulae state that the amplitude of the isodensity statistics depends on the ge- 
ometric and the peculiar velocity effects, whose influences are characterized by the parameter 7. 
The sensitivity to those effects can be understood from the function J, K and L. In figure 2, we 
plot the normalized amplitude 7(7)/ J(l), K{'y,9)/K{l,0) and L{'y,9)/L{l,0) as the functions of 
distortion parameter 7. To show the angle dependence of the functions K and L, we plot the cases 
6 = (solid), tt/6 (long-dashed), it/3 (short-dashed) and it/2 (dotted). Figure 2 suggests that 
the distortion effects may significantly affect the isodensity statistics, especially for the quantities 
described by the function L{'y,9), i.e., G2 and Ni. In next subsection, we will clarify how 7 is 
varied by the distortion effects at a quantitative level. 

Before proceeding the analysis, it should be noted that isodensity statistics in terms of the 
volume fraction / defined by 




instead of the density threshold z^, has been conventionally used in order to compare the Gaussian 
results with non-Gaussian cases (e.g, Gott et al. 1989; Vogely et al. 1994; Canavezes et al. 1998). 
In fact, the isodensity statistics labeled by / remains unchanged under the transformation (35) in 
real space. In cosmological redshift space, however, the statistics cannot be entirely characterized 
by the volume fraction because of the anisotropy. Moreover, we intend to clarify influences of 
the distortion effects including the nonlinear correction. Therefore, we use the straightforward 
expressions labeled by the density threshold v. 
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3.3. Results 

Adopting the above nonlinear model, we now give a prediction for isodensity statistics in 
specific cold dark matter (CDM) cosmological models. In this paper, we consider two typical CDM 
models, in which the cosmological parameters are chosen as {flo, \o,crs,h) = (1.0,0.0,0.6,0.5) 
(SCDM model) and {Qq, Aq, o-g, h) = (0.3, 0.7, 1.0, 0.7) (LCDM model), where erg. the top-hat mass 
fluctuation amplitude at 8/i~^Mpc, is normalized by the cluster abundance (c.g, Kitayama & Suto 
1997). For the unknown biasing factor, we use the time-evolving biasing model by Fry (1996). This 
model simply accounts for the redshift evolution of number density fluctuations in the case that 
the formation and merging process can be negligible. Denoting the amplitude of the biasing factor 
at present by bo, the linearized equation of continuity yields 

K.) = l + ^. (47) 

We fix the value 6o = 1-5 for brevity. 

The distortion parameters 7 and the field correlation (Jq are the most important parameters 
to understand behaviors of the shape of the isodensity statistics. Therefore we first analyze 7 and 
o"o. Figures 3 and 4 show 7 and Uq as the functions of the smoothing radius Rg- Figure 3 shows the 
case adopting the cosmological redshift space with s = SHb(-2) (see eqs.[6] and [8]). The results are 
depicted for the LCDM model ( Upper-panel) and the SCDM model (Lower-panel), though the two 
models show almost similar behaviors. In each model the results at 2: = 1 (Left-panel) and z = 3 
(Right-panel) are shown. 

On each panel the solid line plots the case all the effects are incorporated. Because the solid 
line corresponds to the case observed in the space s, which is referred to 's-space' on the panel. The 
solid line shows that the distortion parameter 7 (the clustering amplitude (Jq) gradually increases 
(decreases) as the smoothing scale Rg increases. To understand the importance of the nonlinear and 

the geometric distortion effects, we also plot the linear results neglecting the nonlinear correction 
(dashed line) and the case in absence of the geometric distortion (dot- dashed line), which are 
referred to 's-space (linear)' and 'q-space' on the panel, respectively. The latter case corresponds 
to the result obtained in the q-space (see eq.[5]), which is the same result obtained in the s-space 
with setting Cy = Cx = 1. 

By comparing the behaviors of the solid and the dot-dashed lines for 7 on large smoothing 
scales, where the linear distortion and the geometric distortion are influential, it is apparent that 
the linear distortion increases 7 and that the geometric distortion decreases 7. On the other hand 
the nonlinear correction becomes important on small scales (compare solid and dashed lines). It 
appears that the nonlinear effects cause more significant influence on the distortion parameter 7 
rather than the clustering amplitude (Tq. This is because the finger-of-God effect compensates the 
gravitational nonlinear growth and the resultant ctq becomes almost same as the linear result. As 
for 7, since the finger-of-God effect dominates the linear velocity distortion on small scales, 7 is 
significantly decreased. It is notable that the correction from nonlinear effects still remains non- 
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negligible even at the larger smoothing length Rg ~ 20/i~ Mpc in the cosmological redshift space 
■Snb- This behavior is caused from the scaling effect of the geometric distortion through the factors 
C|| and Cx- Note also that the difference between the solid line and the dot-dashed line becomes 
significant at higher redshift z = 3, which also indicates the importance of the geometric distortion 
effect at the higher redshift. 

FigTire 4 shows another example, where the cosmological redshift space is chosen as s = SEd(-2) 
(see eqs.[7] and [9]) . In contrast with the case s = Suhiz) in Figure 3a, the geometric distortion 
increases 7 and suppresses the amplitude ctq in the case of the LCDM model ( Upper-panel). For the 
SCDM case, the results in s-space exactly coincides with those in q-space because the cosmological 
redshift space coincide with the comoving real space and the geometric distortion is absent. In 
both cases, correction from the nonlinear effects becomes negligible on the scales Rg ^ 8/i~^Mpc 
at high-redshift. This is also true irrespective of whether the cosmological model of the real space 
is spatially flat or open. 

For comparison, in Figure 3 and 4, the case of no bias is also shown by setting b{z) = 1 {dotted 
line). This alternation does not change the qualitative behavior, however, we can see that 7 and Uo 
are affected by the amplitude of the bias as follows. Increasing the amplitude of the biasing factor, 
7 is decreased, conversely ctq is increased. This effect could be more significant at higher redshift 
because of time-evolution of the bias. 

In principle, behaviors of the isodensity statistics can be entirely understood from Figure 2, 3, 
and 4. As a demonstration, however, in Figure 5 we plot shape-dependence of the genus G3 {upper 
panels) and G2 {lower panels) as functions of the density height u for the LCDM model. Figure 
6 shows the area statistics N^, {upper panels), N2 {middle panels), and Ni {lower panels). In both 
Figures 5 and 6, the cosmological redshift space s = Sin,{z) {s = SEd(^;)) is chosen for the panels in 
the left (right) hand side of the figures. As is expected, the shape of the isodensity statistics is very 
sensitive to the choice of the cosmological redshift space. In the case s = Sui,{z), the significance 
of the nonlinear correction leads to asymmetric shapes of the isodensity statistics even at the large 
smoothing length Rs — 10/i~^Mpc and at high-redshift. The relative amplitude of G2 and A'^i is 
sensitive to the angle 9, which might be a useful tool to determine the cosmological model of our 
universe as a simple geometric test (Alcock Sz Paczynski 1978; see Sec.5). On the other hand, 
in the case s = SEd{z), all the statistics have nearly symmetric shape. This indicates that the 
Gaussian linear results could be a good approximation. Interestingly, at high-redshift z = 3, the 
relative amplitude of the isodensity statistics has the opposite ^-dependence in the case s = SEd{z), 
compared to the one in the case s = SHb(^). As the angle 9 increases, the relative amplitude 
becomes smaller in the case s = SEa{z) and conversely larger in the case s = SHh{z) (c.f. left- and 
right-panels). This behavior comes from the difference of the value 7. For the smoothing length 
adopted here, 7 < 1 in the case s = SHb(-z), while 7 > 1 in the case s = s-Ed{z) (Figures 3 and 4). 
This fact causes the opposite ^-dependence, as is easily deduced from Figure 2. 
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4. Discussion 

The analysis in previous section shows that the anisotropy caused by the geometric and the 
velocity distortions systematically affect the amplitude of the one- and two-dimensional quantities 
with the fractional power by 10 ~ 20%. For the three-dimensional quantities G3 and N^, the 
scaling effect by the geometric distortion enhances or reduces the nonlinearity. These results are 
qualitatively correct as long as a large number of samples are distributed within a narrow range of 
the redshift. In this section, for practical purpose, we focus on the following issues: dctcctability of 
these effects (subsection 4.1); influence of finite redshift range of samples, i.e, the light-cone effect 
(subsection 4.2). 



4.1. Sensitivity 

For quantitative study of high-redshift clustering, we must be careful for the Poisson noise and 
the sampling noise as the main statistical errors. This is not an exception for the case of isodensity 

statistics. While the sampling noise arises from finiteness of the sampling volume, the Poisson noise 
masks the isodensity statistics of underlying density fluctuations when the samples of cosmological 
objects arc sparse and rare. In this subsection, to compare the errors with the systematic effects of 
geometric and redshift-space distortions, we roughly estimate both errors in an analytical manner, 
and investigate sensitivity of the isodensity statistics assuming typical high-redshift samples. 

Let us first estimate the Poisson noise. For simplicity, we consider the isodensity statistics for 
the Gaussian random fields. The analytic formulae in Appendix A reveal that the amplitude of the 
isodensity statistics can be written in a general form : 



(48) 



where i = 3 for G3, i = 2 for G2, and i = 1 for the area statistics, and we defined 



{kl) ^ J— , (49) 

/ d^K W\ksRs)P^"^\ks) 

with the Gaussian window function W{y) = e"^'^/^. In the expression (48), we have simply ignored 
the anisotropy due to the redshift-space distortions. 

The influence of the Poisson noise can be roughly evaluated as follows. The effect of the 
Poisson sampling is incorporated by adding the shot- noise power spectrum (e.g., Feldman, Kaiser 
& Peacock 1994) 

P^^' = (50) 
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where Ug is the mean number density of the observed high-redshift objects. Then we insert the 
following power spectrum P^"'^^{ks) into equation (49), 

p'-'\ks) = p!::^\ks) + p^it' , (51) 

where Ptiuc\ks) means the spectrum free from the shot noise. By substituting equation (51) into 
(48), we can evaluate the ratio of Q^^^ incorporating the Poisson term to the one neglecting it 
(Canavezes et al. 1998). Here we consider the case that the spectrum is given in the power low 
form Ptiu'e' in the cosmological redshift space for simplicity. This yields the correlation 

function of the form 

d^^s D(crd)., X sin(A;,s) 



^ ^ ' J (27r)3 ^ ■ 



kgS V So 



(s) 

where sq is a correlation length. Then the ratio of reduces to 



( 1 + ^ L_ \ 

87r3/2 UsSlRs 

1+ ' ' 



i/2 



1 



V 47r3/2 nsslRs J 



i/2 



where the last line is valid as long as the contribution from the Poisson noise term is small compared 
with the dominant contribution from the Gaussian random field. Thus the quantity defined by 

^"j = S'K^/''nsslRs (53) 

/ shot 

qualitatively gives signal-to-noise ratio resulting from the Poisson noise contribution. 

On the other hand, in estimating the sampling noise, a reliable and robust method is the numer- 
ical technique of bootstrap re-sampling, which has been widely used in the observational/numerical 
study of genus statistics. Although analytical estimation of sampling noise is rather difficult task, 
we can roughly infer the influence of finite sampling effect from number of resolution elements Nj-es 
which provides a useful measure for the statistical power of isodensity statistics (Gott et al. 1989; 
Canavezes et al. 1998). For a given redshift range of the sample region [^min, ^max], we have 

Q 1^ ('^max) S (^^niin)!- 

iVres = ^ , (54) 

where A$7 denotes a fixed solid angle. Gott ct al. (1989) introduced this quantity to discuss the 
quality of data sets and argued that for a fractional accuracy of 25% in measuring the genus curve, 
a sample with A^res = 100 is required. For a measurement with high fractional accuracy of ~ few 
%, large number of resolution element iVres > 10^ is necessary. Note that in a recent study of Point 
Source Catalogue Redshift Survey of IRAS galaxies, number of resolution element reaches more 
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than 300 (Canavezes et al. 1998). In this paper, we use the quantity A^res to estimate the influence 
of the sampling noise. 

Table 1 summarizes the signal-to-shot noise ratio (53) and the number of resolution element 
(54) assuming the typical samples of cosmological objects, SDSS/2dF quasars, Lyman-break galax- 
ies, and cluster samples. Here, we also list the SDSS/2dF galaxies as a representative low-z sample. 
Typical values of the correlation length sq and the mean number density are adopted to evaluate 
(5/A^)shot- To compute the quantity A^res, we assume the cosmological redshift space s — Sb(j(2^) cind 
set the solid angle AO by tt steradian for the SDSS/2dF galaxies, quasars and clusters, 0.25 square 
degree for the Lyman-break galaxies. Prom Table 1, it is clear that the SDSS/2dF galaxy redshift 
surveys unambiguously provide us high quality data preferable to an analysis of the isodensity 
statistics. For our interest of high-redshift samples, the cluster samples can be good objects for 
studying the isodensity statistics. The systematic cosmological effects can dominate the Poisson 
and the sampling noise, and may act as signal. As for the Lyman-break galaxies, while the high 
signal-to-noise ratio {S/N)shot could be attained, the sampling error can dominate the signal due to 
the narrow field of sample region. In an optimistic sense, this is not a problem when the wide-field 
samples are supplied. The photometric SDSS galaxies beyond z ^ 0.2 might be a candidate of 
such samples if the redshifts of objects are precisely measured, although the redshift range will be 
shallower z ^ 1. On the other hand, the quasar samples would be rather sparse and the Poisson 
noise unfortunately dominates the signal, though the sampling error becomes negligible. 

The above estimation of statistical errors might be rather optimistic and idealistic. We do 
not consider the selection effect and the boundary effect of the survey region. Further, we should 
be aware that strength of cosmological effects as well as the statistical errors also depend on the 
biasing of objects, which might be one of the most important obstacle in estimating the sensitivity. 
In a rigorous sense, we cannot give any conclusive statements until elaborating a more quantitative 
study using numerical simulation, however, we can say that the systematic effects of geometric 
and velocity distortions cannot be entirely negligible. The isodensity statistics on clustering of 
high-redshift objects will be measured from the cluster samples, in addition to the high-z galaxies 
from the future wide-field surveys. 

4.2. Light-cone Effect 

Table 1 suggests that observations of a large volume region is more desirable to reduce the 
sampling noise. In this case, the redshift range of the samples becomes large, and the influence 
of the light-cone effect will be important. Namely, the time-evolution of the distribution of the 
objects should be properly taken into account. The light-cone effect on the two-point correlation 
function and power spectrum is rigorously formulated by Yamamoto & Suto (1999) and the effect 
is extensively investigated in various cosmological situations (Nishioka &; Yamamoto 1999, 2000; 
Yamamoto et al. 1999; Suto et al. 2000a, b). 
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Unlike the two-point statistics, which axe computed from the pair-count of objects with a 

fixed separation, the isodensity statistics can be estimated from the pixelated density fields. This 
situation is similar to the analysis of the higher order moments using the count-in-cell method 
(Matsubara et al. 1997). In this CclSG, clS long as the distant-observer approximation is valid, the 
light-cone effect is supposed to be incorporated into the isodensity statistics as follows: 



(-) = """" , ^ , (55) 



^max 




where (i/; z) denotes the i-dimensional genus statistics at redshift z, dV} is the comoving 
volume element in i-dimensional (cosmological redshift) space, and w{z) is a weight function. In 
general, the weight function w{z) arising from the selection effect is determined so as to maximize 
the signal-to-noise ratio. As stated in Suto et al. (1999), the pixelated density field at z is corrected 
by multiplying the inverse of the selection function ^{z) for the density fields and w{z) can be set 
to unity for z^i^ < z < z^^^, in principle. This is actually done when we compute the genus curves 
from the observational data (e.g, Canavezes et al. 1998). Hence, we simply set w{z) = 1. The 
formula (55) can also be applicable to the area statistics by replacing G'f\u; z) with N^''\v; z), the 
z-dimensional area statistics. 

In the above expression (55), the range of the integration [z^i^,z^^^ and the volume factor 
dVc should be determined depending on geometry of the survey depth for the three dimensional 
statistics, and also depending on configuration of a plane (line) underlying the density fields for 
the two (one) dimensional statistics. As an example, let us consider the case that samples are 
distributed in the range of the radial distance s(Zinin) < s < s{z^^) within a fixed solid angle AJ7. 
In this case, the i-dimensional volume factor becomes 

dVP = s{zy-^ds{z) An. (56) 

As a demonstration, in Figure 7, we plot the three-dimensional genus Gg^^' for the Lyman-break 
galaxies (lower panels) and the high-redshift clusters (middle panels) as typical high-z samples, and 
the SDSS/2dF galaxies (upper panels) as a reference. We adopted the parameters of the redshift 
range of the samples and the smoothing radiiis listed in Table 1. The left and right panels show 
the cases s — ^^(^(z) and s — SEd(^)? respectively. For each panel, solid lines show the genus 
curves incorporating the light-cone effect, G'"^'^\i'). On the other hand the dashed and dotted lines 
show the ones on a constant hypersurface at z = z^in and z^ax, i.e., G3 ' (i^; ^min) and G''^\u; z^^^), 
respectively. In the similar way. Figure 8 shows the two-dimensional genus G^^''^ In the case of 
the two-dimensional genus G2, the two-dimensional slice S embedded in a survey field is taken to 
be parallel to the line-of-sight direction, i.e., ^ = 0. 

From Figures 7 and 8, we might regard that the light-cone effect is not significant for these 
samples. That is, as long as the samples listed in Table 1 are concerned, the light-cone effect on the 
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shapes of isodensity statistics is not so significant. This might be because we focused on the shape of 
the isodensity statistics normalized by the ampHtude at the threshold = and/or 1^ = 1, which is 
sensitive to probability distribution function of the random field, rather than that of the clustering 
amplitude. At a high-redshift and/or with a large smoothing radius, the probability distribution as 
a function of the threshold u merely has the weak dependence of time. This consideration suggests 
that the light-cone effect becomes important for smoothing scales on a boundary between linear 
and nonlinear regimes, where the evolution of the probability distribution function can be drastic. 

5. Conclusions 

In the present paper, we have considered various observational effects on the isodensity statis- 
tics of high-redshift objects, and estimated their influences, as a generalization of the work by 
Matsubara (1996). We have presented a simple theoretical model incorporating the geometric dis- 
tortion and the nonlinear correction arising from the gravitational growth and the finger-of-God 
effects, together with useful analytic formulae. Prom this theoretical model (40)-(44), we have 
demonstrated that the geometric distortion affects the shapes of the isodensity statistics and that 
the nonlinear correction can be sensitive to a choice of the redshift-space coordinate as increasing 
the redshift. Based on these results, we have briefly discussed the detectability of these effects by 
evaluating the shot-noisc and the sampling noise contaminations. Furthermore, the influence of the 
light-cone effect has been investigated. The important conclusions are summarized: 

1. When the cosmological redshift space s = Siit,{z) (eqs.[6] and [8]) is chosen, the nonlinear cor- 
rection has an effect on the isodensity statistics for smoothing length larger Rg ~ 10/i~-'^Mpc 
even at the higher redshift z ^ I, where the isodensity statistics have asymmetric shapes. On 
the other hand, when the cosmological redshift space s = s^diz) (eqs.[7] and [9]) is chosen, 
the nonlinear correction turns out to be small and the linear theory of redshift-space distor- 
tions using the Gaussian results (Al), (A2) and (A3)-(A5) could be a good approximation 
for the large smoothing scales. In both cases, the geometric and velocity distortions cause 
the angle-dependence in the low-dimensional quantities, and the relative amplitude can be 
changed by 10 - 20%. 

2. The Poisson sampling error is crucial for the sparse sampling of high-redshift objects like 
the (SDSS/2dF) quasars. The finiteness of the sampling volume might affect the precise 
measurement of isodensity statistics with Lyman-break galaxies. Though more quantitative 
estimation of statistical errors using numerical simulations is needed, however, the isodensity 
statistics can be tested using clusters and high-z galaxies from future wide-field surveys. For 
these samples, the cosmological effects can be comparable or could be dominant to the Poisson 
and the sampling noise. 

3. For a sample that the survey volume is large and the cosmological objects are distributed in 
a wide range of redshift, the light-cone effect can be influential. However, we infer that the 



-19- 



light-cone effect has a weak influence on the isodensity statistics as long as the time evolution 
of the probability distribution function is not drastic, though this picture might depend on 
the evolution of the biasing factor. 

Because the model presented here are quite simple, the predictions should be checked by numerical 
simulations. The prediction can also be sophisticated by using more reliable nonlinear models 
deduced from the simulations. Recent development of cosmological N-body simulation enables us 
to obtain the enormous large volume data over a few Gpc (Colberg et al. 2000). The larger data 
comparable to the Hubble volume size will provide us with chances of quantitative studies for the 
isodensity statistics, as well as the two-point statistics, in which the light-cone effect inevitably 
becomes important. 

Combining our theoretical template with these simulations, various cosmological implications 
might be yielded. One of the most interesting application is the geometric test proposed by Alcock 
&; Paczyhski (1979). As we have seen in Sec. 3. 3, the geometric distortion makes the observed 
structures distorted and it causes deformation of the isodensity contour. Since this effect is sensitive 
to the distance-redshift relation, the anisotropy of isodensity statistics might be a useful tool in 
investigating cosmological models, as well as the anisotropy of two-point statistics (Yamamoto, 
Nishioka &: Taruya 2000). The present theoretical formulae include useful tools for quantifying the 
geometric distortion, as well as for testing statistical nature of biased density fields. To address 
these issues, however, a more quantitative analysis should be developed. 

The uncertainty of the biasing can be a critical issue to clarify the topological features of 
large-scale structure. Throughout the analyses, we have put the linear biasing ansatz (25) and 
the evolution of the biasing simply given by the model (47) neglecting the merging and formation 
processes. In more realistic situation, the biasing of luminous objects is nonlinear and stochastic 
(Dekel &: Lahav 1999) and the merging and formation processes might alter the evolution of biasing 
(e.g, Somerville et al. 2000; Taruya et al. 2000). However, recent theoretical and numerical studies 
suggest that the linear biasing assumption (25) can be a good approximation even in the quasi- 
linear regime (Coles, 1993; Scherrer & Weinberg 1998; Matsubara 1999; Taruya & Suto 2000; 
Taruya et al. 2000). Therefore, in a qualitative sense, wc believe that the present analyses would 
have revealed the generic feature of isodensity statistics on the observed high-redshift objects. As 
for the evolution of biasing, the theoretical prediction based on the Press-Schechter formalism gives 
the biasing factor for the dark halo distribution accurately (Mo & White 1996; Jing 1998; Sheth 
&; Tormen 1999). Incorporating this evolution model into our analytic formulae can provide a 
powerful tool for quantifying the distribution of clusters of galaxies. In any way the isodensity 
statistics may probe a distinctive aspect of the biasing factor. This will be discussed elsewhere 
(Hikage, Taruya & Suto 2000, in preparation). 
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A. Analytic formulae for Gaussian random field 



In this appendix, we give the analytic formulae of isodensity statistics, which have been previ- 
ously derived by Matsubara (1996). Although he only consider the case in absence of the geometric 
distortion (i.e, C|| = c± = 1) by imposing the assumption of linear redshift-space distortion, his 
results can be easily extended to those incorporating the geometric distortion in a rather formal 
way, without any approximation except for the random Gaussianity. 

Let us suppose that the density field is homogeneous random Gaussian. In this case, the 
joint probability distribution function 'P{6^^\ rji, Qj) is divided into 1^(5^^^ Cij)'Pi'ni) regardless of 
the anisotropy. Final expression for the three dimensional genus (30) is immediately obtained by 
following similar calculations in the homogeneous and isotropic case (Bardeen et al. 1986; Hamilton 
et al. 1986). The resultant formula does not alter the shape of the genus as a function of u. Only 
the amplitude is changed^: 

On the other hand, for the two-dimensional genus (31), by using the relations 

[dDiS^""^ - y(yo)8D{rii) |r72sin6' + r/3Cos6'| Cn^ 

= {5d{5^^^ -y(To)C,u) {Soim) |?72sin6l-F?73Cos6'| ), 

1 / \ 2 
-1 /cri,x\ -,,2/2 



and 



one obtains 



(<^d(??i) l^2sin6' + %cos^|) = -^=1(^,9), 



2tt 



Gi'\u,e) = ^ V e-'/^ Lij,e), (A2) 

where the function L is defined by (B4) in Appendix B. In the above expression, we defined the 
distortion parameter 7 given by (see also eq.[45]) 



Similar to the case of the genus statistics, the area statistics, N^, N2, and A^i arc calculated 
using the probability functions being V{6^^\ rji) = 'P{5^^^)V{'qi), which is entirely characterized by 



^This has been partially shown by Matsubara (1996) in the case of the hnear redshift distortion. 
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the field correlations (12) and (13). After some manipulations, one can obtain the final expressions: 



1 <7i,± 
20r (To 



1\pK (Jq 



2i/7r (To 



(A3) 
(A4) 
(A5) 



where L is the same function as defined in (A2) and the functions J and K are also defined in 
Appendix B by (Bl) and (B2), respectively. 



B. Expressions of J(7), ^^(7) and L(7) 

We here present explicit formulae for the functions J, K and L to describe the isodensity 
statistics for Gaussian random field (Appendix A). The functions J, K and L are defined by the 
ensemble average of the Gaussian random field r] as follows: 



V2 



J{l)= — { Jrii + vl + ril ), 
V2 



Hi, 9) = 



K{i, 6)= {yvi + {m sin + % cos O f ), 



{ \r)2 sin 6 + rjs cos 9\ ), 



where the ensemble average (• • •) is evaluated according to the Gaussian probability distribution 
function 



V{r]i)dr]idT]2dT]3 



1 



(27r)3/2(a,.||a2^/2) 



exp 



ah 



2 

m 



driidri2dr]3. 



Prom straightforward calculations, we obtain the following formulae. Note that these formulae are 
also used to obtain the results for the parameterized non-Gaussian case if we replace the quantities 
r]i, cr^ii and afj^ with ([(?3^]^)</, and ( [c?±0]^ )<^, respectively, (see Section 3.2 and Appendix 
C). 



B.l. J(7) 



dxdydz 
(27r)3/2 
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2Q2 

TT 7 



X < 



where we define Q = 7 /(7 — 1 



1 



1 



' Qi^- Q) ~ q7\Q\ 



In 



tan 



1 f 1 



B.2. K(7,^) 



(ItI > 1) 
(ItI < 1) 



(Bl) 



(a;2+j/2+22)/2 



dxdydz 
(27r)3/2 



_ y|B(V(l-72)cos2«) 

where £'(/c) is the complete eUiptical integral of the second kind, 

^7r/2 



/•7r/2 y 

E{k) = d(t> Jl-Psin^,/). 
Jo 



(l7l > 1) 



(l7l < 1) 



, (B2) 



(B3) 



B.3. L{-f,e) 



(27r)3/2 



1 + (7^ - 1) cos 



2 a |V2 



(B4) 



C. Analytic results for the parameterized non-Gaussian case 

In this appendix, we derive analytic expressions for the isodensity statistics with a non- 
Gaussian distribution whose statistical property is characterized by a monotonic function of the 
random field ^(s) (see eq.[35]). Since the field ^(s) is a Gaussian variable with the zero mean and 
the unit variance, the one-point distribution function of the density fluctuations P{S^^^) becomes 



V{6^'^) d5^ 



s) 



1 exp{-i[F-i(<5W)]2} 



2^ 



d5^'^ 
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^ e-1'"/''d<t>. (CI) 



/27r 

Recall that the ensemble average with respect to the density field 5^^^ can be computed through 
the average with respect to the Gaussian field ^. Using the fact that the field variables rji and Qij 
can be written by 

Cij = F,<f> didjcp + F^^ di(j) dj(j), 

the final results are obtained : 

= -I {Sd[F{4>) - va,] 5D[F,^di<t>] 5D[F,^d2<j>] \F,^d^<t>\ 

= -I { ^d[<P - F-\uao)] dD[dicj>] fe[a2</.] \d3^\ {{8idicl)){d2d2^) - {did2cj>f} )^ 

1 fJi 



for the three-dimensional genus (see also Matsubara & Yokoyama 1996), where the ensemble average 
(• • •)<^ is taken with respect to the field (j). For the two-dimensional genus, we have 

G^"V,0) = {5D[F{(t>) - u(To] 5D[F,^d^(t>] \F,^{d24> sm^ + cos^)| 

X [F^ai9i0 + F^<^(ai0)2] >^ 

= -\ {Sd[4> - F-\uao)] 5D[di(f\ |52</> sin^ + 93<^ cos^| Si^k/. ) 



F-i(z.ao)cxp|-i[F-i(i.ao)]'| ^(7), (C3) 



where 7 is given by (45) and we have used the formula (B4) when taking the ensemble average with 
respect to the Gaussian field dicf). As for the area statistics N^, N2 and A^i, repeating the similar 
calculations yields 



' ^ ' ^ 20F ^{F,^m^ 2L ^ / 

iv(^)(j,,0) = 1 ^ '^"^ exp|-i rF-i(z/ao)]'l ^(7). 

^ ^ ' ^ 2Vi ycF>(0F). "^i 2L ^ / 



(C4) 

(C5) 
(C6) 



In deriving the above formulae, we have used the following relations, 

{F{<l>)%, 
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and 

^ = {F4cl>)\{[djct>rU, (C7) 
where the subscript j runs from 1 to 2. 
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Table 1. High-redshift samples and signal-to-noise ratios of isodensity statistics. As a 
representative low-z sample, we also list the SDSS/2dF galaxy samples. Typical values from 
observations are adopted for the correlation length sq and the mean number density Ug. The 

smoothing length Rg is specifically chosen as Rg = 2.5sq except for the quasar samples, 
Rg = 50/i^^Mpc. Then signal-to-noisc ratio (5'/A^)shot is evaluated from definition (53). Number of 
resolution element A^res is computed from (54) assuming the cosmological redshift space s = s^di^) 
and setting the solid angle of survey field AJ7 by tt steradian (SDSS galaxies, quasars and clusters) 

or 0.25 square degree (Lyman-break galaxies). 



samples 




so[^~^Mpc] 


Rs[h~^Mpc] 








SDSS/2dF galaxies 


[0.0, 0.2] 


4.0 


10 


0.1 


710 


3.6 X 10^ 


Clusters 


[0.0, 2.0] 


20.0 


50 


10-5 


9.0 


3.3 X 10^ 


Lyman-break galaxies 


[2.0, 4.0] 


4.0 


10 


10-2 


360 


1.2 X 10^ 


SDSS/2dF quasars 


[0.0, 3.0] 


4.0 


50 


10-6 


0.04 


5.4 X 10^ 
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Observer 



Fig. 1. — A sketch of cosmological redshift coordinate and the configuration of survey region (plane 
S/line L) 
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: K(7.0)/K(1.O) 




7 



Fig. 2. — The normalized amplitude of J {Upper), K [Middle), and L {Lower) as the function 
of the distortion parameter 7 (see also Appendix A). For K and L we fixed ^ = (solid), 7r/6 
(long-dashed), 7r/3 (short-dashed), and 7r/2 (dotted). 
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Fig. 3. — The scale-dependence of distortion parameter 7 and the field correlation ctq as the function 
of the smoothing radius Rs in the case adopting the cosmological redshift space s = SHb(-z) (see 
eqs.[6] and [8]). Upper-panels{Lower-panels) show the LCDM (SCDM) model. The left (right) panel 
shows the case z = 1 {z = 3). Here, we used the biasing parameter given by (47) specifying the 
present value bo = 1.5. The solid and dot-dashed lines indicate the predictions from the nonlinear 
power spectrum (26), where the latter case is obtained by setting Cy = c± = 1. The dashed line 
shows the linear prediction, in which only the linear distortion and the geometric effect arc taken 
into account. For comparison, we plot the result for the no-biasing case, b{z) = 1 (dotted line). 
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Fig. 4. — Same as Figure 3, but we adopt the cosmological redshift space s = s^diz) (see eqs.[7] 
and [9]). Note that the SCDM model (the lower panels) is not affected by the geometric distortion 
because the cosmological redshift-space is exactly same as the real space. Therefore the solid lines 
exactly coincide with the dot-dashed lines in the SCDM model. 
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Fig. 5. — Shapes of the isodensity statistics as the function of the density height v in the LCDM 
model. We assume the cosmological redshift space s = Sm,{z) in left-panel, and s = s^a{z) in 
right-panel. In each panel, we show the result with the redshift and the smoothing scale fixed as 
z = 1 and Rs = 10^~^Mpc (upper-left), z = 1 and Rg = 30/i~^Mpc (upper-right), z = 3 and 
Rs = 10/i~''"Mpc (lower-left), and z = 3 and Rs = 30/i~^Mpc (lower-right). As for the biasing 
parameter, we adopt the Fry's bias model with the present value bo = 1.5 : the three-dimensional 
genus G3 (Upper); the two-dimensional genus G2 (Lower). For G2, we plot the cases 9 = (solid), 
7r/6 (long-dashed), 7r/3 (short-dashed), and 7r/2 (dotted). 
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Fig. 6. — Same as Figure 5, but we plot the mean area of isodensity contour A^3 (Upper), the 
mean length of isodensity contour intersecting with a plane A^2 {Middle), and the mean number of 
crossing points of isodensity contour on a line Ni (Lower). For N2 and iVi, we plot the cases ^ = 0, 
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Fig. 7. — The three-dimensional genus curves G3 evaluated on a light-cone in the LCDM model. For 
each panel, the evolution of bias is assumed to be the Fry's model with 60 = 1.5. The cosmological 
redshift space is chosen as s = Sin-,{z) in left-panel and s = s^diz) in right-panel. The adopted 
range of the redshift-integration and the smoothing scale are shown in each panel. 
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Fig. 8. — Same as Figure 7, but the two-dimensional genus curves G2 evaluated on a light-cone. 
Here, we take the two-dimensional slice S parallel to the line-of-sight direction, i.e., = 0. 



